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Abstract. We extend the shell-model Monte Carlo applications to the rare-earth region 
to include the odd-even nucleus 16I Dy. The projection on an odd number of particles leads 
to a sign problem at low temperatures making it impractical to extract the ground-state 
energy in direct calculations. We use level counting data at low energies and neutron 
resonance data to extract the shell model ground-state energy to good precision. We then 
calculate the level density of I6l Dy and find it in very good agreement with the level 
density extracted from experimental data. 

o ' 

3 ' 1 Introduction 

The shell model Monte Carlo (SMMC) M1I2I3I4I is a powerful method for the calculation of statistical 
properties of nuclei at finite temperature. This approach has proven to be particularly useful in the 
calculation of nuclear level densities B5I6I7I8I91 . Reliable microscopic calculations of the level density 
often require the inclusion of correlations beyond the mean-field approximation. In the SMMC ap- 
proach such correlations are treated in the context of the interacting shell model. SMMC calculations 
can be carried out in model spaces that are many orders of magnitude larger than model spaces that 
can be used in the conventional matrix diagonalization approach to the shell model. 

Recently, the SMMC approach was extended to the rare-earth region and applied to a well-deformed 
even-even nucleus 152 Dy [ 1 1 . Here we extend the application of SMMC in rare-earth nuclei to include 
the odd-even nucleus 161 Dy. The projection on an odd number of particles leads to a sign problem (even 
for good-sign interactions), making it impractical to calculate thermal observables at low temperatures. 
Consequently, the ground-state energy of the nucleus cannot be extracted from direct SMMC calcula- 
tions at very low temperatures. Here we use level counting data at low energies and neutron resonance 
data to determine the ground-state energy. We then calculate the SMMC level density of 161 Dy and 
compare it with the level density extracted from available data. 



2 Choice of Model Space and Interaction 

We use the model space of Ref. iflOl . where the single-particle orbitals are 0^7/2, W5/2, 2si/2, 
O/in/2 and I/7/2 for proton, and O/in/2, 0/29/2, 1/7/2, I/5/2, 2p3/2, 2pi/2, O/13/2, and 1*79/2 for neutrons. 
The single-particle energies are determined to coincide with the spherical Woods-Saxon plus spin-orbit 
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potential in the spherical Hartree-Fock approximation. The effective interaction consists of monopole 
pairing and multipole-multipole terms flOl 

-2 g v P\P v ~YjXa : (O^ + /i; „) • (O x - p + , (1) 

v=p,n A 

where the pair creation operator P\ is given by P\ = 2nZ/m( — y +m+la l, m-v a aj-m-v anc ^ t ' le mu ltipole 
operator <9 i;v is given by 0.i;v = y= Zab(ja\\^Y A \\j h )[al Uv x a ajb - v ]^ with a Jm = {-)j +m a Hn . 
: : denotes normal ordering. The pairing strengths, g v = yg v , where g p = 10.9/Z and g n = 10.9/N 
are parametrized to reproduce the experimental odd-even mass differences for nearby spherical nu- 
clei in the number-projected BCS approximation ItTOl . The quadrupole, octupole and hexadecupole 
interaction terms have strengths given by \a = k^x for A = 2, 3,4 respectively, y is determined self- 
consistently fiTTl and kx are renormalization factors accounting for core polarization effects. Here we 
use the values for y and as in the study of 162 Dy IflOll . 



3 Ground-state Energy 

Because of the sign problem introduced by the projection on an odd number of particles, the thermal 
energy E(J3) of I61 Dy can in practice be calculated only up to/? ~ 5.5 MeV" 1 . The calculations for f3 > 
2 MeV -1 were carried out using the recently implemented stabilization routines [ 10 1 to stabilize the 
canonical propagator. Since the discretization of /3 introduces systematic errors in E(j3), we calculated 
the thermal energy at any given /3 for two values for the time slice (A/3 =1/32 MeV -1 and A/3 — 1 /64 
MeV -1 ) and then performed an extrapolation to A/3 = 0. For /3 < 3 MeV -1 , a linear extrapolation was 
found to be suitable while for larger values of f3, the dependence of E(J3) on A/3 is weaker and we took 
an average value. 

To determine the ground-state energy Eq, we fitted (by a one-parameter fit) the SMMC thermal 
excitation energy E*(T) = E(T) - Eo to the experimental thermal excitation energy. The latter is 
calculated from E*(J3) = -dlnZ(J3)/d/3, where Z is the experimental partition function (see below). 
The experimental partition function and the thermal excitation energy are shown, respectively, by 
solid lines in the top and bottom panels of Fig. [T] The solid circles are the SMMC results shown down 
to a temperature of T a: 0.18 MeV, below which the statistical errors become too large because of the 
sign problem (for comparison the lowest temperature calculated in the even-even nucleus 162 Dy was 
T = 0.05 MeV iflOl ). 

At sufficiently low temperatures the experimental partition function can be calculated from Z(J3) = 
2,(27, + l)e~P E < where Ej are the experimentally known energy levels of the nucleus (measured with 
respect to the ground-state energy). This partition function and the corresponding thermal excitation 
energy are shown by the dashed lines in Fig. Q] However, since the level counting data is incomplete 
above a certain excitation energy, the average experimental thermal energy becomes saturated at tem- 
peratures above T ~ 0.18 MeV (see the inset in the lower panel of Fig . fT} . Since SMMC data exist only 
above T ~ 0.18 MeV, a realistic estimate of the experimental thermal energy at higher temperatures is 
necessary to determine Eq from the fit. We accomplish this by using an experimental partition function 
defined by 

N ^oo 

Z(T) = YpJi + Ve- E,/T + I dE PBBF (E)e- E/T , (2) 

,■ ^ En 

where En is an energy below which an essentially complete set of levels is known and pbbf(E) = 
j£a~ l t 4 (E - Ay 5 ^e 2 V fl (-E-^) j s an experimentally determined level density parametrized by the back- 
shifted Bethe formula (BBF) 1 12]. The single-particle level density parameter a and the backshift A 
are determined from level counting data at low energies and the s-wave neutron resonance data at the 
neutron resonance energy. We find a = 18.564 MeV -1 andzl = -0.615 MeV. The solid lines in Fig.[T] 
are then calculated from Eq. (f2). The number of the low-lying states (determining En) is chosen 
such that the two curves (solid and dashed) for the experimental partition function merge smoothly at 
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Fig. 1. Partition function Z (top panel) and thermal excitation energy E* (bottom panel) versus temperature T 
for 161 Dy. The SMMC results (solid circles) are compared with the results deduced from experimentally known 
levels (dashed lines), and from Eq. (|2} (solid lines). The ground-state energy Eo is obtained by a one-parameter 
fit of the SMMC thermal energy to the experimental thermal energy (solid curve). Inset: thermal energy E versus 
temperature for T < 0.4 MeV. 



a sufficiently low temperature. We found that N = 8 is a reasonable choice. Fitting the solid curve in 
the top panel of Fig.[T]to the SMMC thermal excitation energy (solid circles), we found a ground-state 
energy of Eo = -363.920 + 0.020 MeV. We emphasize that the SMMC partition function and thermal 
energy fit very well the corresponding experimental quantities by adjusting only a single parameter 
/••(I- 



4 Level Density 

The level density is the inverse Laplace transform of the partition function. Its average is determined in 
the saddle-point approximation in terms of the canonical entropy and the heat capacity, which in turn 
are calculated from the average thermal energy E(J3) [5|. The SMMC level density is shown by the 
solid circles in Fig. [2] We compare it with the experimentally determined BBF level density (dashed 
line). The level counting data at low energies are shown by the histograms, and the neutron resonance 
data is shown by the triangle at E x = 6.454 MeV. We find that at low excitation energies and close to 
the neutron resonance energy, the SMMC results agree quite well with the experimental level density. 
At intermediate energies, the SMMC level density is slightly below the BBF level density. 



5 Conclusions 

We have extended the SMMC calculations in the rare-earth region to include the odd-even nucleus, 
161 Dy. The sign problem introduced by the projection on the odd number of particles makes it im- 
practical to calculate directly the ground-state energy. We have circumvented this problem by using 
the neutron resonance data and level counting data at low energies. We have then calculated the level 
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Fig. 2. Level density of 161 Dy. The SMMC results (solid circles) are compared with experimental results. The 
histograms are from level counting data, the triangle is the neutron resonance data and the dashed line is the 
experimentally extracted BBF level density. Inset: level density for excitation energies E x < 1.5 MeV. 



density and find it to be in good agreement with the experimental level density extracted from avail- 
able data. The method outlined here requires sufficient level counting data at low energies and neutron 
resonance data. 
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